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Abstract: We study noncommutative field theories at finite temperature to learn 
more about the degrees of freedom in the non-planar sector of these systems. We 
find evidence for winding states. At temperatures for which the thermal wavelength 
is smaller than the noncommutativity scale, there is a drastic reduction of the de- 
grees of freedom in the non-planar sector. In this regime, the non-planar sector has 
thermodynamics resembling that of a 1 + 1 dimensional field theory. 
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1. Introduction 

Gauge theories on noncommutative spaces can be obtained by taking the infinite tension 
limit of string theory in the presence of a very strong B-field [1, 3, 4, 2]. One can thus 
look upon these field theories as interesting limits of string theory in such an extreme 
environment. Regardless of their ancestry, these field theories have very interesting 
properties [5]-[28], making their study fascinating in its own right. Many unusual 
properties already emerge in scalar field theories. In particular, one sees the emergence 
of unexpected infrared behavior in correlation functions [4] (we will refer to this paper in 
the following as MSV) unlike anything seen in conventional field theories. These long 
distance correlations appear in massive theories and are not related to any massless 
fields in the Lagrangian. Rather, these infrared effects are due to the Moyal phases. As 
can be seen in Feynman diagrams, it is the integration over the high ultraviolet modes 
that is responsible for this peculiar infrared behavior. The UV and IR are intertwined 
in a manner never seen in local field theory on commutative spaces. MSV have a 
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proposal for what kind of degrees of freedom are responsible for these IR correlations: 
their claim is that in the zero slope limit and large B field, the closed string states do 
not decouple. 

To gain some more insight into the structure of these noncommutative field theories, 
we subject them to a heat bath and consider various thermodynamic quantities. We will 
show that the partition functions of these theories get contributions at high temperature 
from configurations winding around the temporal direction. This implies that there are 
extended objects in these theories that have the ability to wind, which seems consistent 
with the suggestions of MSV. We also test this picture by compactifying the spatial 
commuting direction to a circle and again find contributions to the partition function, 
this time from configurations winding around the spatial circle. 

It is useful to consider a simple field theory like massive 4 and the supersymmetric 
Wess-Zumino model. We were motivated to include a supersymmetric example in order 
to show that the existence of winding states persists in theories where the sensitivity to 
the UV physics is softer. Because of the special UV properties of the supersymmetric 
case, the IR sensitivity is milder and allows us to limit ourselves to the lowest order in 
perturbation theory when calculating thermodynamic quantities. 

The paper is organized as follows: in section 2 we will describe in some detail the 
thermodynamics of noncommutative 4 . There we present the calculation of the free 
energy to 0(g 2 ) and discuss the various divergences and also show how the winding 
states appear. In section 3, we investigate whether supersymmetric systems whose UV 
behavior are less singular affect the existence of winding states. For this purpose, we 
consider the Wess-Zumino model and find again the existence of winding states. We 
calculate the free energy, from which the entropy, internal energy and specific heat can 
be derived. In section 4, we return to the 4 theory and show how to calculate the free 
energy in the "mean field" approximation. In section 5, we discuss the implications of 
these results and briefly repeat the calculations at zero and finite temperature for the 
case where the sole commuting direction is a circle. We show that at finite temperature, 
there are no membrane states wrapping around the torus spanned by the temporal and 
x 3 cycles. Instead, we find winding states winding each of the cycles independently. 

We then conclude with some speculations. 

2. The thermodynamics of noncommutative g 2 4 in perturba- 
tion theory 

We present a perturbative calculation of the free energy, F, which assumes that the cou- 
pling constant g 2 is small. As we will see, the appearance of IR divergences invalidates 
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the perturbative expansion. Nevertheless, we believe it is instructive to go through 
this exercise as it already shows the existence of winding states. Later in section 4 we 
will have to use a mean field approximation to ensure that long range correlations are 
properly screened. 

In this paper, we will restrict ourselves to the case where the noncommutativ- 
ity is solely among the spatial directions.. Without loss of generality then in this 
4-dimensional case, we can choose the (1,2) plane to be noncommutative and leave 
time and the third spatial dimension commutative. 

At finite temperature, the Feynman rules are the ones that correspond to a theory 
with "time" being a circle of radius (3, where (3 is the inverse temperature, (3 = 1/T. 
This implies that the integrals over frequencies in diagrams are replaced by discrete 
sums over the so-called Matsubara frequencies. The noncommutativity in the (1,2) 
plane manifests itself in Feynman diagrams through the appearance of Moyal phases 
in vertices. These Moyal phases depend on momenta in the (1,2) directions. 

When evaluating the contributions to the free energy, one needs to consider both 
planar and non-planar diagrams. The planar contributions are very much like in the 
commutative theory except for some combinatoric factors [5] . Inspite of the changed 
combinatorics, the behavior in the UV is still controlled by zero temperature physics. 
In order to illustrate the divergences that occur in the perturbative expansion of the 
free energy in the noncommutative theory, we will first concentrate on the non-planar 
contributions to lowest order in g 2 . This will also enable us to show the appearance of 
contributions from winding states to the free energy. 

The leading non-planar contribution at order g 2 to the free energy comes from a 
two-loop diagram. 



Fig.l: Planar and nonplanar two loop contributions to the free energy. 

This diagram is divergent, where some of the divergences can be interpreted as 
originating in the UV and some divergence can be attributed to the IR. The UV diver- 
gences will be shown to arise from zero temperature physics and similar UV divergences 
will appear in higher order terms in the expansion in g 2 . Also, in a manner similar to 
the IR divergences that appears at lowest order, each successive order of perturbation 
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theory will be plagued by IR divergences. These IR divergences grow more singular 
with increasing order in the perturbative expansion, invalidating the expansion. We 
will return to the formal application of this mean field approximation in section 4 after 
the study of the 0(g 2 ) contribution to the free energy in the Wess-Zumino model. 
The contribution to the free energy from this two-loop non-planar diagram is: 



00 00 r d 3 r> rl 3 h p ip9k 

-a 2 T 2 V V / P — — - (2 1) 

9 J^oo (2?r) 3 (27T)3 (W +p2 + M2)( 4^ + F + M2) > 



where pdk = 9(p\k2 — p 2 ^i) and p 2 = p 2 + p 2 2 + p^ 2 . 

We will be interested in temperatures much larger than the mass, so that to leading 
order in f3M, it can be neglected. 

The integration over the momenta can be done by introducing Schwinger parame- 
ters. In order to show how contributions of winding configurations appear, we will look 
at equation (2.1) in some detail. The integral appearing in this expression can after 
some manipulations be rewritten as: 

_„i T Y I -^L - (2 2) 

ZjJ (27r) 3 (W +p2)(47r 2( Z2/ 32 + |0p|2) +A -2)' ^ • > 

where \9p\ = 9(pi 2 + p 2 ) 1 ^ 2 and A is an ultraviolet cutoff similar to the one used by 
MSV. We also have neglected the mass in obtaining this expression since we are working 
in the regime f3M <C 1. 

The sums over n and / can be performed and will lead to the following formula 
involving Bose distributions 



2 f d 3 p 1 + 2^1)1 + 2^(27^1) 
9 J (2vr) 3 2\p\ 4n\6p\ ' 1 ' 

where \p\ = {p\ 2 +P2 2 +P^ 2 ) 1 ^ 2 and np(\p\) = g \ l \ _ x is a Bose distribution at temperature 
T = 1/(3. 

Note the appearance in the last expression of two Bose distributions, one at tem- 
perature T and the other at temperature 1/T. This would be a rather bizarre state 
of affair for a system in thermal equilibrium. One possibility is that the system is not 
in equilibrium, but why then in this very special way, with one set of the degrees of 
freedom distributed thermally at temperature T and the other at temperature 1/T? 
We will choose the option that the system is indeed in equilibrium, and that there are 
configurations, if one thinks in a path integral approach to the partition function, that 
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wind around the temporal circle. If one looks at the integral appearing in equation 
(2.3), one can see the contributions of winding states with momentum \9p\ to the free 
energy. 

Another interesting feature of equation (2.3) is that the p 3 = sector is invariant 
under the interchange of /3 and 2^9/(3. This transformation looks like T-duality where 
the string scale l s 2 is replaced by 9. We should remind the reader though that this 
property was obtained in the f3M <C 1 limit. We do not know whether this duality 
invariance persists beyond this order to the full theory. 

The expression in equation (2.3) has UV divergences and IR divergences. The 
IR divergence is linear and comes from the small (^1,^2) region of integration. This 
divergence, as will be shown later, is cured by the mean field approximation. 

The UV divergent terms in the non-planar sector are calculated in the appendix. 
Apart from a zero temperature contribution to the vacuum energy, there is a linearly 
divergent term which is also linear in the temperature. This is a rather unusual contri- 
bution to the free energy. Indeed, this term does not contribute to the internal energy 
of the system but does provide a temperature independent contribution to the entropy. 
This contribution to the entropy is due to the degeneracy of the winding states. As 
can be seen in the Bose distribution, their dispersion relation does not involve the mo- 
mentum along the commuting spatial dimension, p 3 . The very large number of winding 
states with momenta p± t 2 = are the source of the linear divergence. This divergence 
disappears when screening of the large distance correlations is taken into account. The 
super symmetric case does not have such a divergence at the leading order in g 2 because 
of its milder behavior in the IR. 

3. The thermodynamics of the noncommutative Wess-Zumino 
model 

The Lagrangian for the noncommutative Wess-Zumino model is: 

£ = ii) llV( T>' v + A*D A - -Mipip - -M^ - gip* r.l - #* c .1 F*F , 
where F is given by 

F = -MA* - gA* * A* . 

The calculation of the free energy to 0(g 2 ) proceeds along the same lines as in the 
bosonic case. In this case there are additional contributions due to new interactions 
among the bosons and bosons with fermions. The fermions, as usual, have frequencies 
that are odd multiples of the temperature. 
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Three pairs of diagrams contribute at order g 2 , each pair consisting of a planar and 
a non-planar piece. They sum to give 



2 f d 3 p d 3 k l + e ipdk f . , , ,x , , , . ,x 
F/V = -g J — — ——{n B (u p )+n F (u p )) (n B (u k ) + n F (u> k )) 

M 2 

+ 9 2 T 4 0(— ) (3.1) 

with cj p = Vf 2 + ^ 2 and tib,f(u p ) = l/(e f3u}p =F 1). Notice that this expression for the 
free energy at high temperature (f3M <C 1) is "T- duality invariant" as in the bosonic 
case. This property applies to the sector with zero momentum along the commuting 
direction. This result can be seen by expanding one of the Bose distributions in a series 
of exponentials, 

-| oo 

/3n(p 2 +M 2 ) 1 / 2 



1 71=1 



e 



g/3(p 2 +M2) 1 /2 

Then performing the integral over d 3 p, neglecting the mass, produces the series, 

1 



E 



n ^n 2 [3 2 + \9k\ 2 ' 

A similar expression can be obtained for the integration over the Fermi distributions, 
except that in the fermion case the series is an alternating one. The presence of winding 
states is then clear as well as the "T- duality invariance" of the free energy. 



Fig. 2: Planar and non-planar two loop fermionic contributions to the free energy. 



Fig. 3: Planar and non-planar two loop bosonic contributions to the free energy. 
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The expression (3.1) for the non-planar contribution to the free energy can be 
further evaluated to yield 



' dr > L d 1 „:„u„ ( 3 - 2 ) 



y/2 0(3 2 Jo Jo sinh q 

At high temperature, (3M <C 1, we consider two limiting cases, T 2 <C 1/6 1 and 
T 2 ^> 1/0. 1 The first limit corresponds to the case where the thermal wavelength is 
bigger than the length scale associated to the noncommutativity. In this limit, one 
might expect to recover conventional high temperature behavior. Indeed, the leading 
non-planar contribution to the free energy here is 

vl~-» 2T '- < 3 - 3 > 

In the other limit, T 2 ^> 1/9, the thermal wavelength is smaller than the noncom- 
mutativity scale. We therefore expect novel behavior. We find the free energy in this 
regime to be 



F 
V 



,T 



2 



»„ -!>T logTe ' (3 ' 4) 

where the logarithm was obtained by numerical methods. In this case, the contribution 
to the free energy resembles the contribution of a massless gas in 1+1 dimension, 
modified by a logarithm. 

One is tempted to conclude from these results that there is a drastic reduction of 
degrees of freedom in the non-planar sector! For thermal wavelengths larger than the 
noncommutativity scale, the system behaves as a relativitic 3 + 1 dimensional gas. For 
wavelengths smaller than this scale, the system behaves roughly like a 1 + 1 dimensional 
gas. As the temperature is raised past 1/9 1 / 2 it looks like the winding states do behave 
as strings. Up to the logarithm, one is reminded by the temperature dependence in 
equation (3.4) of the result obtained by Atick and Witten [29] for the high temperature 
behavior of strings. We do not know whether this is just a coincidence or whether 
there is some deeper significance to this behavior. A test of this high temperature 
dependence of the free energy due to the non planar sector would be to study higher 
dimensional examples where there are more noncommuting spatial dimensions. It is 
amusing to note that the non-planar sector of this system exhibits, to leading order 
at very high temperature, the equation of state p = p, p being pressure and p energy 
density. 



1 We also assume that 9M 2 <C 1, i.e. the Compton wavelength of the scalars is larger than the 
length scale associated to the noncommutativity. 
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4. The thermodynamics of noncommutative g 2 (j) 4 in the "mean 
field" approximation 

We now go back to the noncommutative g 2 <p A theory of section 2. There, the appear- 
ance of IR divergences ruins the validity of the perturbative expansion (in g 2 ) for the 
free energy. This is a familiar phenomenon in thermal physics. For example, in the 
case of a plasma of electric charges in QED, one finds IR divergences in perturba- 
tion theory for the free energy. This can be traced to the large distance behavior of 
the correlation function, (p(x)p(0)). This IR divergence disappears after summing a 
series of diagrams contributing to the correlation function between charge densities, 
(p(x)p(0)). This series of diagrams is a geometric series in the one-particle irreducible 
charge density correlation function. When this series is summed, one finds that the 
charge densities are screened. The screening mass or Debye mass is the IPI charge 
density correlation function evaluated at zero momentum. 

In the case of noncommutative 4 , the IR divergence already appears at zero tem- 
perature and is again screened by summing a geometric series in the one-particle irre- 
ducible two-point function (MSV). 

In the finite temperature case, we will approach the problem of IR divergences by 
treating the gas in a mean field approximation. In this approximation, each quantum 
of the scalar field is moving in the average field produced by the other scalars. Unlike 
more traditional cases where the mean field approximation is used, this case has the 
particular feature that the mean field will depend on the momentum of the particle 
that moves in this background. This momentum lies in the noncommutative plane. 
Using this approximation, we will show how the long range correlations are screened. 

The Hamiltonian H is: 

H=^Jd 3 x (n 2 (x) + (V0(a:)) 2 + M 2 <p 2 (x)) + g 2 J d 3 x </>(x) * <j>(x) * <p{x) * <j>(x). 

The free energy F is obtained from the logarithm of the partition function Z, where 
Z = tre _/3// , as — f3F = logZ. When g 2 is small, it can be formally written as Fq, the 
free gas contribution to F, corrected by 0(g 2 ) contributions. We should remind the 
reader that this procedure fails because of IR divergences. We present this expression 
for F here because it will be useful when we get to the mean field approximation: 

V" = / Wf [ 2 + M2)V2 + l0g(1 " ex P<-^ 2 + M2 ) V2 ] + °^)- 

The contribution from the free gas to the previous formula could have been obtained 
by calculating a one-loop diagram in the path integral formulation. 
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In the mean field approximation, a picture emerges where the system looks like a 
free gas for which the dispersion formula for the energy of a quantum is modified by 
the average effect of the other quanta. 

Formally, this will be implemented in the path integral by using the following 
action: 

d 3 k 



S 



dr 



(d t <f>Y(k) - k l 4>\k) - M 2 <f) 2 (k) - g 2 (P) ek p (k) 



(4.1) 



where (<t> 2 )ek satisfies an integral equation: 

d 3 p 



(0 2 V = 



(2vT 



,3 4?r 2 n 2 



+ p 2 + M 2 +g 2^ ep 



We will not solve this integral equation exactly but rather find the solution to it as a 
power series expansion in g 2 . In this paper we will confine ourselves to work with the 
solution to leading order. 

In the limit of high temperature, M/3 < 1, we neglect to first order the mass M 
and we sum over n, and obtain to leading order (4> 2 )ok- 



<0V 



E 



i 



^ 4ir 2 n 2 [3 2 + (2ir\6k\) 2 + A- 2 

' (1 + 2n 1 / j3 (J4ir 2 \6k\ 2 + A~ 2 ) 



2/3^(2vr|^|) 2 + A" 



(4.2) 
(4.3) 



Notice that the Bose distribution in the expression for ((f) 2 ) is unusual as it dis- 
tributes k\ t 2 according to a "temperature" r, where r = 1/9T, i.e. proportional to 
the inverse temperature. One can interpret this property, as was argued earlier, to 
be evidence for winding states. Indeed, the summation that appears in (0 2 ) is over 
an integer multiple of the length of the temporal circle rather than the inverse radius 
as would be the case for momenta. The free energy can then be written as a "tree" 
contribution plus a one-loop contribution. 



d 3 p 



PL = M 2 (j d A X(f) 2 (x)) + / d 4 xq)(A ) *o(.r) * o(.r) * o(.r)) + 

2 



lp 2 + M 2 + g 2 (<P 2 ) 9p + log (l - exp{-(3{p 2 + M 2 + g 2 (<p 2 )e P f 2 
In this last equation, the average value of 

M 2 J d 4 x (j) 2 (x) + ^- J d 4 x <f>(x) * (f)(x) * <j>{x) * <j>{x) 



(4.4) 
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is evaluated by using the action S, equation (4.1), in a path integral calculation. 

As expected and demonstrated explicitly in the appendix, the IR divergences dis- 
appear. 

The one-loop contribution 



y p 2 + M 2 +g 2 {(j) 2 )ep + lQg _ exp ( _ /9(p 2 + M 2 + g^^) 0p f^ 

can be visualized as a single quantum of the scalar field propagating around a loop with 
a momentum that at 0(g°) is thermally distributed at temperature T. This particle 
is moving in the presence of a background, (0 2 ), which it interacts with at 0(g 2 ) and 
that consists of states winding around the temporal circle of radius 1/T. These winding 
states have momenta solely in the noncommutative plane. It is quite tempting to think 
of them as closed strings reminiscent of the suggestions made in MSV. 

To be complete, one should study the UV divergences and see whether they are 
associated to zero temperature effects. We will not attempt this study here. It should 
be kept in mind that it is not completely clear yet whether this theory is renormalizable, 
even though some evidence exists. 




5. Winding states around spatial cycles 

The presence of winding states, which in the finite temperature case wind around 
the temporal circle, can also be noticed when one compactifies the commuting spatial 
dimension, x 3 . The states winding around rr 3 can be detected for example if we calculate 
the free energy for this case. 

The calculation of the previous section can be repeated, with each / ^ replaced 
by where 1 is an integer. 

To show the existence of winding states, it is sufficient to again consider the integral 
equation for the two point function: 

( , 2) _L^f^p ^ 

W)9k PLjjJ (2vr)2^ + 4^ +p2 + M2 + ^ (02)ep - 

In the limit of high temperature, f3M <C 1, and small radius, LM <C 1, we find for 
(4> 2 )ek, upon performing the / d 2 p, 



L2 



)ok - 



4ir 2 n 2 (3 2 + Att 2 PL 2 + (2ir\6k\) 2 + A~ 
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This last expression clearly shows the presence of winding states around both cir- 
cles. Again these states become prevalent in the small radii limit, which is to be 
expected since in these limits these states are very "light" . 

This formula for {4> 2 )ek also suggests that there are no wrapping states around the 
2-torus spanning the temporal and x% dimensions, which would look like an additional 
piece (nl f3 L) 2 in the denominator of this equation. 

The same result holds in the supersymmetric case. In some ways, this could have 
been expected. Imagine subjecting this system to a thermal environment. In this case, 
the bosons and fermions have different distributions and we should not therefore expect 
a drastically different behavior from the non-supersymmetric cases. 

One might speculate about the existence of multidimensional extended objects in 
higher dimensional theories where several % ^ 0. To test this, we briefly consider the 
example of the six- dimensional <f> 3 theory, with various 9^ 7^ 0. 

At high temperature and after compactifying the commuting x^, on a small circle, 
it should be possible to see whether or not wrapping states emerge. This can be tested, 
as we learned above, simply by calculating ((p 2 )e P ,e' P '- We obtain 



We thus conclude that there are no wrapping states, just winding states. 
6. Conclusions 

What is the physics that is responsible for the existence of these winding states? 

Clearly the UV-IR connection that exists in the noncommutative plane is crucial. 
The expressions which reveal the existence of the winding states show that this peculiar 
IR behavior originates in the UV region of integration. One might then conclude that 
the uncertainty relationship between x\ and x 2 in the case where 9 12 7^ is implemented 
by or implies the existence of one-dimensional structures. In other words, measuring 
x\ sharply renders the other direction x 2 totally uncertain and maybe this is intimately 
related to the existence of extended one-dimensional objects in a rather mysterious 
way. These "strings" can then wind around various cycles. However, we did not find 
multidimensional objects wrapping multicycles in higher dimension. 

We have discovered from studying the high temperature regime of the supersym- 
metric theory the existence of two different behaviors. We believe that these results are 
valid in general, but we have only been able to calculate reliably to leading order in the 
supersymmetric case. When the thermal wavelength is bigger than the noncommuta- 
tivity length scale, we find that the non-planar sector behaves as a conventional 3 + 1 




l.n 
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dimensional relativistic gas. For the thermal wavelength smaller than the noncommu- 
tativity scale, we observe a reduction of the degrees of freedom in the non-planar sector. 
The non-planar sector now behaves as a 1 + 1 (not, as might have been expected, 2 + 1) 
dimensional gas. This seems to indicate a drastic reduction of the degrees of freedom in 
that sector. Equivalently, we find the equation of state of this sector is approximately 
that of an incompressible gas with equation of state p = p. 

It would be interesting to see whether this behavior persists in higher dimensional 
theories. 

In the theories we have studied, the fascinating behavior of the non-planar sector 
is overshadowed by the planar sector. Thus, the challenging question remains of what 
theory captures exclusively the physics of the non-planar sector. 

Acknowledgments 

The work of WF, EG, JG, AK-P, SP, PP is supported in part by the Robert Welch 
Foundation and the NSF under grant number PHY-9511632, JG is partially supported 
by AEN 98-0431 (CICYT), GC 1998SGR (CIRIT), SP is also supported by NSF grant 
PHY-9973543. 

A. UV divergences in g 2 (j) 4 

Of the four terms in equation (2.3), the two damped by the Bose distribution ra^d^l) are 
obviously UV convergent. The remaining two give rise to a linear and a logarithmic UV 
divergence (we remind the reader that the argument \9p\ of the second Bose distribution 
nyp appearing in equation (2.3) is independent of ps). 
The term 

_ 2 T f d3 P ni/p(\0p\) 
9 J (2tt) 3 \p\ 4ir\6p\ 

is divergent in the region p 3 ^> \jp\ + p\. In this region, we obtain 

Notice that the x integral is IR divergent. If we introduce an IR cutoff to regulate 
this divergence, the substitution we use to extract the temperature dependence of the 
expression renders this cutoff temperature dependent. We ignore this problem at this 
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point, as we expect it to be an IR artifact that will be cured upon resummation of the 
appropriate diagrams. 

The linear divergence arises as: 



d 3 p 1 1 g 2 T_ A 



(2vr) 3 2|p|47r|^| 4(2tt) 2 9 

As discussed above, this divergence contributes to the entropy of the system, but not to 
the internal energy. We therefore attribute it to the degeneracy of the winding states 
as discussed in the text. 



B. IR divergences in g 2 cj) A 

The only potential infrared divergence in the mean field expression (4) for the free 
energy stems from the term 

f3V( J d 4 x(f)(x) * (f)(x) * <f)(x) * <p(x)). 

The non-planar contribution to this term is obtained from equation (2.1) by replacing 
the bare by the dressed propagators: 



_ 2 T 2 V f d 3 p d 3 k 

9 h J W ( 2n ) 3 + ^ + MH 5 2 (^} flj ,)(^ + F + MH ' 

(B.l) 

where (0 2 ) is given to leading order by equation (4.3). For small p and small k, both 
terms in the denominator of the above expression blow up and no infrared divergence 
can occur. For small p and large k (or vice versa), we can ignore (<p 2 )ek and perform 
the d 3 k integral as in (2.2). Upon performing the sums, we obtain for small p 



2 f _fp_ 1 + 2MH) 2 
J (2?0 3 2y/p* + M* + g*(<f?)e p ** 2 \°P\ 2 ' ' 

As one factor of \0p\~ 1 is absorbed by the (<p 2 )e P under the square root, this expression 
is finite in the infrared. 
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